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<N ■ ABSTRACT 

^ ■ We follow the nuclear reactions that occur in the accretion disks of stellar 

mass black holes that are accreting at a very high rate, 0.01 to 1 Mq sec , as 
is realized in many current models for gamma-ray bursts (GRBs). The degree 
^ ■ of neutronization in the disk is a sensitive function of the accretion rate, black 

^ ■ hole mass, Kerr parameter, and disk viscosity. For high accretion rates and 

Tj- . low viscosity, material arriving at the black hole will consist predominantly of 

neutrons. This degree of neutronization will have important implications for 



O 



CN . the dynamics of the GRB producing jet and perhaps for the synthesis of the r- 

^ . process. For lower accretion rates and high viscosity, as might be appropriate for 

Q_|. the outer disk in the coUapsar model, neutron-proton equality persists allowing 

Q . the possible synthesis of ^^Ni in the disk wind. ^^Ni must be present to make 

any optically bright Type I supernova, and in particular those associated with 
S : GRBs. 
> '■ 

^ . Subject headings: gamma rays: bursts — nucleosynthesis — accretion disks 

1. INTRODUCTION 



Growing evidence connects GRBs to the the birth of hyper-accreting black holes (Fryer, 
Woosley, Sz Hartmann 1999), that is stellar mass black holes accreting matter from a disk 
at rates from ~0.01 to 10 MQsec~^. Such models include the collapsar (Woosley 1993; 
MacFadyen Sz Woosley 1999), merging neutron stars and black holes (Eichler et al. 1989; 
Janka et al. 1999, 2001), supranovae (Vietri & Stella 1998, 1999), merging hehum cores 
and black holes (Zhang & Fryer 2001), and merging white dwarfs and black holes (Fryer et 
al. 1999). For such high accretion rates the disk is optically thick, except to neutrinos, and 
very hot, consisting in its inner regions of a (viscous) mixture of neutrons and protons. At 
its inner boundary the disk connects to the black hole. Whatever processes accelerate the 
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putative GRB-producing jet might therefore be expected to act upon some mixture of disk 
material and other background medium (e.g., the collapsing star in the collapsar model). 
Farther out, material will be lost from the disk in a vigorous wind (MacFadycn & Woosley 
1999; Narayan, Piran, & Kumar 2002). In fact, Narayan et al. suggest that most of the 
disk will be lost to a wind except for those models and in those regions where neutrino losses 
dominate the energy budget. 

It is thus of some consequence to know the composition of such disks. In the least 
case, the nucleosynthesis may be novel and could account for rare species in nature, such 
as the r-proccss. At most, the presence of free neutrons may affect the dynamics of the 
GRB jet (Dcrishev, Kocharovsky, & Kocharovsky 1999; Fuller, Pruet, & Abazajian 2000), 
the GRB neutrino signature (Bahcall & Meszaros 2000), hght curve (Pruet & Dalai 2002), 
and afterglow (e.g. via a "pre-acceleration" mechanism similar to the one discussed by 
Beloborodov 2002). It is also of some consequence to know whether the disk wind consists 
of radioactive ^^Ni, as is necessary if a visible supernova is to accompany the GRB. If the 
electron mole number. Ye — T,{ZiXi/Ai), is less than 0.485, the iron group will be dominated 
by ^^'^^Fe and other more neutron-rich species (Hartmann, Woosley, & El Eid 1985) which 
will be incapable of illuminating the supernova. Since the jet itself is inefficient at heating 
sufficient matter to temperatures required for nuclear statistical equilibrium (T z5 x 10^ K), 
supernovae seen in conjunction with GRBs (Galama et al. 1998; Bloom et al. 2002, for 
example) would be difficult to understand. 

We have thus undertaken a survey of the nucleosynthesis that happens in the disks of 
rapidly accreting black holes. The work is greatly facilitated by the existence of numerical 
and semi-analytic solutions that yield the temperature-density structure and drift velocity 
(Popham, Woosley, & Fryer 1999). These solutions have been verified in the case of the 
collapsar model by direct numerical simulation (MacFadyen & Woosley 1999). 



2. COMPUTATIONAL APPROACH 

For a disk that is optically thin to neutrinos the neutron to proton ratio, n/p, is deter- 
mined by the competition between electron and positron capture on nuclei and free nucleons. 
Since the baryon number per co-moving volume is conserved, it is convenient to work with 
the electron fraction Ye defined by 

Ye^E Zi{Xi/Ai) (1) 

where Zi, Xi, and Ai are the proton number, mass fraction, and atomic mass number (integer) 
of the species i. We shall refer to compositions with Ye < 0.5 as "neutron-rich". 
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If we neglect lepton capture on bound nuclei and approximate the disk material as 
consisting of a mixture of free nucleons and a-particles, the evolution of with radius is 
given by 



u . (VFe) = V{r)^ T^JVWcT^ ^ "^^-^ r \ ' ^ J ■ 

(2) 

Here X„ is the mass fraction of free nucleons, u is the 4-velocity of the flow, M represents 
the black hole mass, and V{r) is the radial drift velocity as measured in an inertial frame co- 
rotating with the disk. In the middle term in Eq. 2 we have assumed steady state conditions 
as well as cylindrical symmetry, and we have adopted the relatively simple Schwarzschild 
metric relevant for most of our calculations. Generally, relativistic effects do not make a 
big difference in our calculations of the electron fraction, though they are important for 
deter ming the structure of the disk. 

In Eq. 2 Ae-p and Ae+n are the rates for the processes 

e~ + p ^ n + Ue (3) 
e+ + n p + Pe- (4) 

These rates are given by 

Ae-p = K w\w-—fG.{l,w)S4T,UF)dw (5) 
Jsrn rrie 

me 

Ae+„ = K w^{w + — fG+{l,w)SJTMdw. (6) 

Here K 6.414 • 10~^s~^ determines the free neutron lifetime, 5m ~ 1.293 MeV is the 
neutron proton mass difference, and m^ is the electron mass. The functions and 5'+ 
are the electron and positron distribution functions, while the G± are the Coulomb wave 
correction factors discussed in Fuller, Fowler, & Newman (1982). 

The electron Fermi energy C/p is found by inverting the expression for the net electron 
number density, 

n: -nt^ pYeNA = ^ i^ff [ P\S-{T, C/p) - 5+(T, U^))dp (7) 
with p — s/w'^ — 1. 

We solve for the evolution of Y^. by integrating eq. 2 for the results presented in Popham, 
Woosley, & Fryer (1999). Those authors assumed Ye — 1/2. Because we will show that in 
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some cases Yg <S 1/2, and consequently that the electron Fermi energy and degeneracy pres- 
sure are substantially smaller than for = 1/2, our results are not entirely self consistent, 
but should suffice given the very approximate nature of a one-dimensional calculation. 

3. RESULTS 

3.1. Scaling with Disk Viscosity and Accretion Rate 

We present results for several different disk models. The parameters describing these 
models, as well as the electron fraction near the event horizon at r ^ lO^cm, arc given in 
Table 1. In all cases we assume Ye = 1/2 at large radii. For simplicity we have adopted 
a constant black hole mass of 3 Mq though the scaling relations for mass are obvious in 
Popham et al. All else being equal, a larger mass will give less electron capture. 

Table 1 shows a clear trend with disk viscosity and accretion rate: electron capture 
in the inner disk becomes more pronounced with decreasing viscosity and increasing mass 
accretion rate. This arises because low viscosity flows inefficiently advect angular momentum 
outwards and are, for a given mass accretion rate, denser than high viscosity flows. Dense 
flows are electron degenerate, so that pair creation is suppressed and electron capture 
dominates over positron capture. Similarly a larger accretion rate also implies a denser disk 
- for a given value of a - and more capture. 

Before discussing in some detail the evolution of Yg in our calculations we note that the 
influence of black hole spin on the composition is also important. This is because the angular 
momemtum imparted to the black hole via the accreting matter will typically drive the Kerr 
parameter high. For example, a Kerr parameter of a ~ 0.9 is typical for collapsars. A larger 
Kerr parameter allows the disk to move to smaller radii before entering the event horizon, so 
that more electron capture and a smaller Yg will result. This is demonstrated by model G, 
which corresponds to a black hole with a = 0.95. In this model n/p near the event horizon 
is about ten times larger than n/p for the same model except with a — (model B). The 
influence of the black hole spin on Y^ is limited to regions quite near the event horizon. For 
r > 10^-^cm, Ye in model B is essentially the same as in model G. These considerations imply 
that as far as the composition is concerned, the main difference between flows around zero 
angular momentum holes and those around large angular momentum holes will be that the 
jet in the large angular momentum case will likely be more neutron rich. The composition 
in the bulk of the wind coming from the disk, however, will likely be similar in the two cases. 

In Fig. 1 we show the evolution of Yg with radius for an accretion disk with a = 0.1 
and a relatively low mass accretion rate (at least for standard coUapsar models), M — 
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0.01 Mq sec~^. This flow is not dense enough to drive the electrons degenerate. Consequently 
Yg is governed by thermal capture. In this case the only asymmetry is the neutron-proton 
mass difference and actually increases slightly owing to the threshold for the rate in Eq. 3, 
becoming greater than 0.5. The free nucleon mass fraction appearing in Fig. 1 is calculated 
from the simple estimate provided in Qian & Woosley (1996). Lepton capture on heavy 
nuclei is negligible for this disk. 

Results for a hotter and denser flow, M = 0.1 MqScc"^, a = 0.1, are shown in Fig. 5. 
The material at radii r < 10 ''•^cm is mildly degenerate and the electron fraction is driven 
to Ye ~ 0.44. The quantity \e-p{r/V) is the product of the electron capture rate and a 
rough measure of the dynamic timescale, or time left before the infalling material crosses 
the event horizon. Because Y^ can not come to equilibrium unless \e-p{r/V) ~ 1, Y^. is not 
in equilibrium in this flow. Also, because Xe-p{r/V) <^ 1 when bound nuclei are present, 
weak processes on heavy nuclei are not important even in the extreme limit where the bound 
nucleons behave as free nucleons with respect to electron capture. 

Figures 3 and 4 show the influence of viscosity on the composition of the flow. For 
a — 0.03 electron capture becomes important at r ~ 10^'^cm and Yg is in close equilibrium 
with capture until r ^ 10^'^cm. There is at most ~ 1/10 of an electron capture per bound 
nucleus in the disk. For a = 0.01 the inflowing material becomes degenerate and neutron rich 
early on. Most of the free protons are locked into a-particles until the a-particles dissociate. 
Positron capture on the excess free neutrons results in a brief increase in Y^., as seen in the 
bump at r ~ lO'^'^cm in Fig. 2. The final n/p in this disk is very large, ~ 20. 

The influence of the mass accretion rate on the composition is seen in Fig. 5, where 
we plot results for M = 1 M0 sec^^ and a = 0.1. This disk is quite similar to the M = 
0.1 MqSCC"^, a = 0.3 case. Indeed, this is just what is expected from the scaling relations 
in Popham, Woosley, & Fryer (1999). For M ^ 1 Mq sec^^ and larger, the assumption that 
the disk is optically thin to neutrinos begins to break down (MacFadyen & Woosley 1999; 
Di Matteo , Perna, & Narayan 2002). This influences both the structure of the disk, and, 
because neutrino capture becomes important, the evolution of Yg. Because neutrino trapping 
likely influences the structure of the disk at the same level as the assumption that Y^ — 1/2, 
we do not calculate for disks calculated under the assumption of partial neutrino trapping 
(and Ye = 1/2). 
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3.2. Effects of Neutrino Capture 

The effect of neutrino capture on tlie evolution of can be addressed in an approximate 
way. In particular, we are interested in the rate Ai/^n for the process 

i/e + n^p + e". (8) 

The rate for capture is suppressed both because of the low proton number density and the 
absence of high energy positrons in the flow. The number of neutrinos captured per neutron 
is roughly T^e,cap"'i/e/"'n ~ ^e7"i.e,cap"'i.e/"'p> where T^^,cap is the optical depth for the process 
in Eq. 8, and nv^jn^ ('^i^e/'^p) is the number of electron neutrinos produced per neutron 
(proton) in the flow. The ratio R of the number of neutrino captures per neutron to the 
number of electron captures per proton is 

R^^= T,,,cap>;. (9) 
'^e-p 

This equation implies that an optical depth of 1 to capture is roughly equivalent to a 
doubling of Ag+n- We note here that while the the fe's produced in electron-capture reactions 
are not thermal, their average capture cross section is only approximately a factor of 2 
higher than the average capture cross section calculated under the assumption of a thermal 
distribution for the electron neutrinos (at the electron Fermi energy and temperature) for 
conditions of interest in accretion disks. 

It is difficult to quantitatively calculate the influence of neutrino captiu^c on Y^. in a 
post-processing step. This is both because of the complexities of neutrino transport, as well 
as because of the feedback between Fe, the electron capture rate, and the disk dynamics. 
However, the following simple considerations argue that neutrino capture is unlikely to have 
a dramatic influence on Y^. 

To discuss some specific cases, consider first the M = O.lMQsec"^, a = 0.1 disk. Di 
Matteo , Perna, & Narayan (2002) argue that neutrinos will have an absorptive optical 
depth of unity at r ~ 3 ■ lO^cm for this disk. Because weak processes have essentially frozen 
out by this radius (see the Ae-p(r/y) curve in Fig. 2), neutrino capture will have a small 
influence on Yg. 

Denser flows trap neutrinos more efficiently but arc still expected to remain neutron rich. 
For example, for the M = IMqScc^^, a = 0.1 disk, Di Matteo , Perna, & Narayan (2002) 
show that the neutrino absorptive optical depth will be greater than unity for r < 2 • lO^cm. 
As weak processes are rapid compared to the dynamic timescale in this disk, the result 
will be an increase in Y^. To estimate the magnitude of this increase we adopt the simple 
procedure of artificially increasing the positron capture rate, Ae+n — 5Ae+n, which is likely 
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an overestimate of the protonization rate. The increase in Ag+nresults in an electron fraction 
of Yf. PS 0.166 near the event horizon. This is approximately 50% higher than the electron 
fraction calculated by neglecting the influence of i/^ capture, though still quite neutron rich. 

To summarize, Ye in disks with high viscosity and modest accretion rates will not be 
affected by neutrino capture simply because neutrinos are not trapped or because weak 
processes are too slow. The electron fraction in disks with neutrino absorptive optical depths 
of a few will increase somewhat but will remain neutron rich owing to the sharp rise in Ag-p 
with Yg. 



3.3. Transport to the Surface of the Disk 

In order for a low electron fraction to have observable implications the electron fraction 

must not come to equilibrium at Y^, = 1 /2 as the material travels out of the plane of the 
disk. To estimate the evolution of the electron fraction in convective blobs moving out of the 
disk we parametrize the convective timescale by the turbulent convection speed (Narayan, 
Piran, & Kumar 2002) 

Vtmh ~ OiVK — 6.3 X lO^a-irf^^^ cm (10) 

where q;_i = a/ 0.1 and ry = r/lO^cm. This implies a time to go one pressure scale height 

Tconv ~ r/vtnrh ~ 16 aZl ms. (11) 

Wc also assume that the convective blob expands adiabatically. An estimate for the entropy 
per baryon in the disk is given by Qian & Woosley (1996): 

s/h ~ 0.052:^ + 7.4 + In ^ . (12) 

pio y Pio ) 

Here Tmbv is the temperature in MeV and pio is the density in units of 10^°gcm~^. The first 
term on the right hand side of Eq. 12 is the contribution to the entropy from relativistic hght 
particles (7/e^), while the next two represent the contribution from free nucleons. Eq. 12 is 
not appropriate when electrons are degenerate. For degenerate electrons the contribution of 
thermal pairs is small and the coefficient 0.052 should be closer to 0.02, representing just 
the entropy of the photons. However, changing the coefficient to 0.02 makes little difference 
for our purposes. 

When relativistic particles dominate the entropy, the adiabat satisfies oc p. When 
free nucleons dominate the entropy, the adiabat satisfies T^/^ oc p. The evolution of the 
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electron fraction in a convective blob is crucially sensitive to the adiabat on which the blob 
travels. 

To illustrate this, consider first a disk with M — 0.01, a — 0.1. In Fig. 6 we show 
the evolution of the electron fraction in matter originating from two different points in this 
disk and for different assumptions about Tconv The upper curves in Fig. 6 correspond to 
material originating from r = 10^'^cm in the disk. At this point the entropy is relatively 
high, s/kh ~ 23 by the estimate in Eq. 12, and consequently the adiabat is closer to oc p 
than to T-'^/^ oc p. The temperature decreases quite slowly relative to the density, allowing 
for the formation of pairs which efficiently drive n/p to equality. For large Tconv this process 
occurs more efficiently than for small Tconv The lower curves in Fig. 6 correspond to matter 
originating from r = lO^cm in the disk. At this point the entropy is ~ 17, so the adiabat is 
near T^/^ oc p. For low entropies, then, the temperature and density decrease at a comparable 
rate and pair formation is somewhat suppressed, keeping Yg low. 

For flows where the electron fraction is high enough to drive the inner disk neutron 
rich, the entropy is generally dominated by non-relativistic particles, and the adiabat again 
preserves a low Yg. This is illustrated in Fig. 7 where we plot the evolution of electron 
fraction in an adiabatic convective bubble for two different disks characterized by a low 
Yf,. In this figure we only show the evolution of a fluid element from a single initial radius 
(10 ''cm) because the entropy in these is disks never dominated by relativistic particles when 
Yg is low, so that the behavior shown in Fig. 7 is generic. 

The above considerations outline the general features of how Y^ evolves in material as 
it travels from the center to the surface of the disk. In flows where Y^ is driven small, the 
degeneracy of the electrons implies that free nucleons dominate the entropy. A low Ye then, 
will remain low. By contrast, in flows where Y^ is close to 1/2, the adiabat can be close 
enough to T"^ oc p for Yg to be driven to 1/2, and the neutron excess at the center of the disk 
is larger than the neutron excess in material flnding its way out of the disk. 

4. IMPLICATIONS 

4.1. Radioactivity in the disk wind 

The most obvious consequence of our calculations is that ^^Ni will be absent from the 
winds of hyper-accreting black holes unless a) the accretion rate is low ( sO.l Mq sec~^), and 
b) the disk viscosity is high, a ^0.1. Interestingly, modern views regarding a-disks and GRB 
models favor values close to these. 
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Typical accretion rates in the coUapsar model are 0.05 - 0.1 MqSCC"^ (MacFadyen & 
Wooslcy 1999, their figures 5 and 10). The neutron excess will be smaller in Type II coUapsars 
powered by fall back rather than direct black hole formation (MacFadyen, Woosley, & Heger 
2001). Considerable accretion into the hole and mass loss from the disk may continue, 
at a declining rate, even after the main GRB producing event (~20 s) is over (Zhang & 
Woosley 2002). It thus seems likely that the coUapsar model will be able to provide the 
^^Ni necessary to make the supernovae that accompany GRBs (though only if a ^0.1). This 
is also true of the slower accreting models like helium-star black hole mergers and black hole 
white dwarf mergers. However any wind from merging compact objects, or similar models 
like the supranova, will be neutron rich. Though perhaps of interest for nucleosynthesis, 
they will not produce ^^Ni, at least during the black hole accretion epoch. 



4.2. Neutron excess in GRB jets 

For disks with high accretion rates, and certainly for merging neutron stars or black 
hole neutron star mergers, the matter near the event horizon will be very neutron rich. If 
this material pollutes the outgoing jet, the the GRB jet will itself, at least initially, contain 
free neutrons. 

The dynamics of accelerating neutron-rich jets can differ dramatically from the dynamics 
of pure proton jets. Fuller, Pruet, & Abazajian (2000) showed that a high Lorentz factor 
fireball that is neutron rich can lead to two very distinct kinematic components, a slow 
neutron outfiow and a fast proton outfiow. This arises because the uncharged neutrons 
are weakly coupled to the radiation dominated plasma, and are accelerated principally via 
strong neutron-proton scatterings (Derishev, Kocharovsky, & Kocharovsky 1999). Strong 
scatterings freeze out when they become slow compared to the dynamic timescale and at 
this point the neutrons coast. If this decoupling occurs while the jet is still accelerating, 
then the coulomb-coupled protons go on to have a larger Lorentz factor than the neutrons. 

Roughly, dynamic neutron decoupling is only expected for fast jets. Here the precise 
meaning of "fast" depends, among other things, the fireball source size. For relativistic fiows 
originating from compact objects, the dynamic timescale characterizing the acceleration of 
the flow is about 1 ms and the final Lorentz factor must be greater than ~ 300 in order 
for dynamic neutron decoupling to occur. For jets in the coUapsar model, the timescale 
characterizing the acceleration is set by the surface of last interaction of the jet with the 
stellar envelope at ~ lO^^cm. In this case the final Lorentz factor of the jet has to be > 3000 
for dynamic neutron decoupling to occur. Because this is an impossibly high Lorentz factor 
for the coUapsar model, dynamic neutron decoupling will not occur. 
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There may be a detectable neutrino signature of dynamic neutron decoupling. Neutron- 
proton collisions occurring during decoupling will generate pions. In turn these pions will 
decay and lead to the generation of neutrinos with energies of a few GeV in our reference 
frame. These neutrinos should be detectable at the rate of about one per year in next 
generation neutrino telescopes (Bahcall & Meszaros 2000). In addition, there may also be a 
direct electromagnetic signature of neutron decoupling. For some bursts arising from external 
shocks, a slow decoupled neutron shell will decay and shock with the outer proton shell as 
the outer shell plows into the interstellar medium. This leads to a characteristic two-peaked 
structure in the burst (Pruet & Dalai 2002). 

Regardless of whether dynamic neutron decoupling occurs, neutrons in the jet may have 
observable implications for the GRB afterglow. This is particularly true when the shocking 
radius, i.e. the radius at which the shocks giving rise to the observed 7— rays occur, is smaller 
than the length over which free neutrons decay. This condition is 

^shock < 72TnC ~ 10^^72Cm, (13) 

where rghock is the shocking radius, 72 is the Lorentz factor of the outflow in units of 100, 
and Tn ~ 1000 sec is the free neutron lifetime. Eq. 13 is generally satisfied for bursts from 
internal shocks and bursts from external shocks in a very dense medium. When Eq. 13 holds, 
the neutrons can stream ahead of the slowing proton shell (the complement of the dynamic 
neutron decoupling discussed above) and deposit energy as they decay. Implications for the 
resulting afterglow in such a case have been discussed by Beloborodov (2002). 

In the next section we discuss the r-process in winds and jets from accretion disks. How- 
ever, we first note that GRB jets are characterized by interesting light element synthesis. 
Pruet, Guiles, & Fuller (2002) and Lemoine (2002) calculated the thermal synthesis of light 
elements in GRB-like outfiows. They find that the final deuterium mass fraction can be of 
order 1%, some three orders of magnitude larger than the primordial yield. For kinemati- 
cally well coupled flows, the nucleosynthesis depends sensitively on Y^., with deuterium yield 
decreasing with increasing neutron excess. For flows in which dynamic neutron decoupling 
occurs, high energy neutron-a collisions will spall deuterons and can result in flnal deuteron 
mass fractions as high as 10% (Pruet, Guiles, & Fuller 2002). 

4.3. The r-process 

There are two possible sites for the r-process here in those cases where Yg is low - in the 
jet and in the disk wind. 
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4.3.1. In the jet 

The jet has the merit of originating in the vicinity of the black hole where the neutron 
excess is likely to be greatest. If magnetic fields drive the heating and initial acceleration 
of the jet, then this neutron excess will likely be preserved. If neutrinos drive the outfiow 
(MacFadyen & Woosley 1999), Yg will be reset to some extent by the weak interactions. 
However, if the details of the outflow above the black hole are similar to spherically symmetric 
neutrino driven ultra-relativistic winds, then only i/P annihilation will be important and the 
flow will remain neutron rich (Pruet, Fuller, & Cardall 2001). 

It is clear that the rapid expansion of the jet will be favorable to freezing out with a 
large abundance of free neutrons and this is conducive to the r-process (Hoffman, Woosley, 
& Qian 1997). However, a pure nucleonic jet coming from the adiabatic expansion of a 
fireball with energy loading 7] = Eiatemai/ pc^ ^100 is too much of a good thing when it comes 
to entropy and rapid expansion. The entropy was given in eq. 12 (Qian & Woosley 1996) 

V^ = ^(^)Wp (14) 



The internal energy (erg/gm) is 



60 ^ he 

and hence, for r) = 200 and kT ~ 2 MeV 

W^g-^j^ ^ ^ 

The time scale for the expansion, r/c, is less than a millisecond. The high entropy 
implies a low density at the temperature when n and p can start to recombine. Coupled 
with the rapid expansion, one might conclude that the jet will remain pure nucleons and 
light nuclei and that an r-process is impossible. 

This conclusion would be erroneous on several counts. First, the jet may not always 
have such high energy loading. There may be cases where no gamma-ray burst is produced 
and the energy loading of the jet is much less. The entropy would then decrease and the 
expansion time scale would increase. Whether to call these "dirty fireballs" or simply an 
extension of the "disk wind" is simply a matter of taste (see below) . 

Of more relevance to GRB models, and particularly to those that involve massive stars, 
is that the jet will not escape the star without interaction. Along the walls of the jet, the 
Kelvin Helmholtz instability will occur (Zhang & Woosley 2002). The jet will also encounter 
one or more shocks where it will be abruptly slowed. Most interesting is the "jet head" which 
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moves through the star subrelativistically with a speed that varies but is of order c/3. At 
this head the jet mixes with stellar material and the mixture is swept backwards (in the 
moving frame) forming a cocoon. It typically takes 5 to 10 sec for the jet head to penetrate 
the roughly one solar radius of overlying Wolf-Rayet star in the coUapsar model. 

During this time, neutrons in the jet are mixed with heavy nuclei - He, O, Si, Ne, Mg 
and the like - which can serve as seeds for an r-process. The overall dynamics is likely to be 
quite complex and its study would requires a multi-dimensional relativistic simulation well 
beyond the scope of this paper. However, Zhang & Woosley (2002) find typical densities 
and time scales in the jet prior to break out arc 10^"* to 10^^ g cm~^ and seconds. This is 
enough that many, if not all of the neutrons in the jet would capture. 

Still, the overall yield in gamma-ray bursts is too small to account for the solar r- 
process. Assume there is one GRB for each 100 supernovac. Each burst has 10^^ - 10^^ 
erg of relativistic ejecta (both jets) and a Lorentz factor that is at least 200 (Lithwick & 
Sari 2001). The jet will be some mixture of entrained stellar material and jet, but assume 
that the initial jet is half or more of the relativistic ejecta. This gives an equivalent yield of 
about 10~^ Mq per supernova. Probably the r-process will be at most 10% of this (the rest 
may be a-particles and entrained nuclei). This then is several orders of magnitude less than 
required to make the solar r-proccss. However, GRB jets might still be important sources of 
the r-process in metal-deficient stars or of select rare nuclei in the sun. 



4-3.2. The r-process in the disk wind 

A more significant r-process could come from the non-relativistic ejecta that comprise 
the disk wind, provided that the entropy is increased by magnetic processes to well above its 
value in the middle of the disk, as in for example, Daigne & Mochkovitch (2002). Assume 
that the accretion is 50% efficient, that is that half the material that enters the disk eventually 
enters the black hole and the other half is blown away (Narayan, Piran, & Kumar 2002). 
A typical collapsar powered GRB involves the accretion of from one to several solar masses, 
so ~1 Mq is ejected from the disk. Most of this will not be neutron- rich enough or expand 
fast enough to make the r-process, but consider the implications if only 1% of the mass that 
is lost comes from the inner disk at times when the accretion rate is over 0.1 M0sec"^ and 

n < 0.4. 

So close to the hole, significant entropy will be added by any acceleration process that 
produces a strong outfiow. A crude estimate is s/A;;, ~ GMm„/ {rkTmax) ~ 400 r^^(2 MeV /Tmax) 
(Qian & Woosley 1996), where T^nax is the temperature where the energy deposition is the 
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greatest. In fact the conditions in the inner disk for accretion rates ~0.1 MqSCC"^ are not 
so different - in terms of neutrino luminosity, neutron excess, temperature, and gravitational 
potential - from the neutrino-driven wind in an ordinary supernova, and one might expect 
a similar r-process (Woosley et al. 1994). The gravitational potential though, can in prin- 
ciple, be greater and a larger entropy and shorter time scale might be realized. Assume 
that the ejecta will be very rich in a-particles and contain of order 10% r-process by mass. 
Putting the numbers together then gives an r-process synthesis equivalent to 10~^ Mq per 
supernova, making this possibility well worth further investigation. 
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Table 1. Electron Mole Number Near the Event Horizon 



Model 






Ye 


A 


0.01 


0.1 


0.510 


B 


0.1 


0.1 


0.435 


C 


0.1 


0.03 


0.119 


D 


0.1 


0.01 


0.045 


E 


1.0 


0.1 


0.115 




0.03 


0.1 


0.527 




0.1 


0.1 


0.077 



^Accretion rate in Mq sec ^ 
*^Disk viscosity 

^Disk properties inferred from 
the scaling relations in Popham, 
Woosley, & Fryer (1999) for this 
model. 

•^Kerr parameter a = 0.95 for this 
model. 
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Fig. 2. — Evolution of Yg and X„ in a disk with M = O.IM© sec~^ and a = 0.1. Shown also 
is the quantity Xe-p{r/V), a rough measure of whether or not there is time for to come 
to equilibrium with respect to capture. 
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Fig. 5. — Same as Fig. 2, except for M — IMq sec ^ and a — 0.1. 
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Fig. 6. — Evolution of electron fraction in a convective and adiabatic fluid element in a disk 
with M — O.lM0sec~^ and a — 0.1. The assumed convective timescale Tconv and radius in 
the disk from which the material originates is given next to each curve. 
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Fig. 7. — Evolution of electron fraction in a convective and adiabatic fluid element in two 
different disks characterized by low electron fractions. The assumed convective timescale 
and disk parameters are given next to each curve. In all cases the material is assumed to 
originate from r = lO^cm. 



